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Abstract 

Charged jet production cross sections in p-Pb collisions at = 5.02 TeV measured with the 
ALICE detector at the LHC are presented. Using the anti-A:T algorithm, jets have been reconstructed 
in the central rapidity region from charged particles with resolution parameters R = 0.2 and R = 0.4. 
The reconstructed jets have been corrected for detector effects and the underlying event background. 
To calculate the nuclear modification factor, Rppb, of charged jets in p-Pb collisions, a pp reference 
was constructed by scaling previously measured charged jet spectra at ^ = 7 TeV. In the transverse 
momentum range 20 < px.chjet < 120 GeV/c, /?ppb is found to be consistent with unity, indicating 
the absence of strong nuclear matter effects on jet production. Major modifications to the radial 
jet structure are probed via the ratio of jet production cross sections reconstructed with the two 
different resolution parameters. This ratio is found to be similar to the measurement in pp collisions 
at = 7 TeV and to the expectations from PYTHIA pp simulations and NLO pQCD calculations 
at ^iNN = 5.02 TeV. 


©2015 CERN for the benefit of the AEICE Collaboration. 

Reproduction of this article or parts of it is allowed as specified in fhe CC-BY-4.0 license. 


See Appendix lAl for the list of collaboration members 






Charged jets in p-Pb 


ALICE Collaboration 


1 Introduction 


Jets are the observable final state of a fragmenting parton produced e.g. in scattering of partons in 
nuclei with a large momentum transfer, Q^. At sufficiently large Q^, the jet production cross section 
is computable since it can be factorized into the non-perturbative parton distribution and fragmentation 
functions and the cross section of partonic scatterings, which is calculable in perturbative QCD (pQCD) 
lU]. Jet measurements in p-Pb and their comparison to pp provide a tool to better constrain effects of 
(cold) nuclear matter on these factors. In particular, they can be used to examine the role of a modification 
of the initial distribution of quarks and gluons, e.g. shadowing effects and gluon saturation Jz yfl, and 
the impact of multiple scatterings and hadronic re-interactions in the initial and final state 10, Ho . 

In central heavy-ion collisions, the production of jets and high-pj particles is strongly modified: in Pb-Pb 
collisions at the LHC, the observed hadron yields are suppressed by up to a factor of seven compared to 
pp collisions, approaching a factor of two suppression at high Ig-M]. A similar suppression is also 
observed for reconstructed jets in central Pb-Pb |l9-13]. This phenomenon, referred to as jet quenching, 
has also been observed previously in high-px particle production in central Au-Au collisions at RHIC 
UMi. It is attributed to the creation of a quark-gluon plasma (QGP) in the final state, where hard 
scattered partons radiate gluons in strong interactions with the medium as first predicted in 21]. This 
results in an radiative energy loss of the leading parton and a modified fragmentation pattern. 

Initially, p-Pb collisions have been seen as the testing ground for isolated cold nuclear matter effects, 
without the formation of a hot and dense medium. However, recent results on low-px particle production 
and long range correlations in p-Pb collisions at = 5.02 TeV 122142511 exhibit features of collective 
behavior, similar to those found in Pb-Pb collisions, where they are attributed to the creation of a QGP. 
At high px, results on the production of unidentified charged particles 1126142911 and jets 11301 l3in in p-Pb 
collisions at = 5.02 TeV are consistent with the absence of a strong final state suppression. The 
question to what extent other nuclear effects lead to an enhancement of particle production at high px is 
still open, a possible enhancement in p-Pb collisions has been reported for single charged hadrons 0281] . 
The measurement of jets in p-Pb collisions compared to single hadrons tests the parton fragmentation 
beyond the leading particle with the inclusion of low-px and large-angle fragments. 

A jet is defined experimentally by the algorithm that combines the measured detector information such 
as tracks and/or calorimeter cells into jet objects and by the parameters of the algorithm. The desired 
properties of such algorithms in pp(p) collisions and in the corresponding theoretical framework have 
been discussed e.g. in 021] . In general, jet algorithms aim to reconstruct the kinematic properties of the 
initial parton with as little dependence on the details of its fragmentation process as possible, i.e. the al¬ 
gorithms should yield consistent results when applied in a theoretical calculation at any stage of a parton 
shower and at final state particle level. A particularly well suited class of algorithms in this context are 
those using sequential recombination schemes, which are infrared and collinear safe, in contrast to many 
conceptually simpler cone algorithms. The computationally optimized implementation of sequential re¬ 
combination algorithms in the FastJet package 0311 facilitates their applicability also in collision systems 
with high multiplicity and thereby the comparison of results obtained with the same jet algorithms in 
pp, p-Pb, and Pb-Pb collisions. An additional complication in the context of jet reconstruction in high- 
multiplicity events arises from the large background particle density, i.e. particles in the same aperture 
as the jet that are not related to the initial hard scattering. This background can be subtracted on an 
event-by-event basis and the impact on the reconstructed jet observable needs to be evaluated carefully 

iShS- 


In this paper, jets reconstructed from charged particles {charged jets) with the anti-kx algorithm measured 
with the ALICE detector in p-Pb collisions at -^^nn = 5.02 TeV are reported for different resolution 
parameters, R. Section 2 describes in detail the correction steps needed in the analysis, including the 
effect of the event background and its fluctuations on the jet observables and the unfolding procedure to 


2 























Charged jets in p-Pb 


ALICE Collaboration 


account for background as well as detector effects. The results are presented and discussed in Sec. 3. 


2 Data Analysis 

2.1 Event and Track Selection 


The data used for this analysis were taken with the ALICE detector Il36l] during the p-Pb run of the 
LHC at = 5.02 TeV at the beginning of 2013. Minimum bias events have been selected requiring 
at least one hit in both of the scintillator trigger detectors (VOA and VOC) covering the pseudorapidity 
2.8 < Tjiab < 5.1 and —3.7 < T7iab < —1.7, respectively 0711 . Here and in the following, T]iab denotes the 
pseudorapidity in the ALICE laboratory frame. Compared to this frame (with positive rj in the direction 
of the VOA), the nucleon-nucleon center-of-mass system moves in rapidity by jnn = —0.465 in the 
direction of the proton beam OSH . 


The event sample used in the analyses presented in this manuscript was collected exclusively for the beam 
configuration where the proton travels from VOA to VOC (clockwise). A van-der-Meer scan was used 
to measure the visible cross section Ovo = 2.09 3=0.07 b for this case 13911 . Monte Carlo studies show 
that the sample consists mainly of non-single diffractive (NSD) interactions and a neg ligible contribution 
from single diffractive and electromagnetic interactions (for more details see 1 3^, ^]). The trigger is not 
fully efficient for NSD events. This inefficiency affecfs only evenfs wifhouf a reconsfrucfed verfex, i.e. 
wifh no particles reconsfrucfed wifhin fhe accepfance of fhe SPD. The loss of efficiency is esfimafed fo 
be 2.3% wifh a large sysfemafic uncerfainfy of 3.1% ll3^ . In fhis paper, fhe normalization to NSD events 
is only used for the construction of the nuclear modification factor. 

In addition to the trigger selection, timing and vertex-quality cuts are used to suppress pile-up and bad 
quality events. The analysis requires a reconstructed vertex, which is the case for 98.2% of the events 
selected by the trigger. In addition, events with a reconstructed vertex |z| > 10 cm along the beam axis 
are rejected. In total, about 96M events are used for the analysis. 


Charged particles are reconstructed as tracks in the Inner Tracking System (ITS) 114Ih and the Time Pro¬ 
jection Chamber (TPC) which cover the full azimuth and |T]iab| <0.9 ll42ll . Lor tracks with reconstructed 
track points close to the vertex (from the two inner Silicon Pixel Detector (SPD) layers of the ITS), a 
momentum resolution of 0.8% (3.8%) for pT = I GeV/c (50 GeV/c) is reached ll^ . The azimuthal 
distribution of these high quality tracks is not completely uniform due to inefficient regions in the SPD. 
This can be compensated by considering in addition tracks without reconstructed track points in the SPD. 
Lor those tracks, the primary vertex is used as an additional constraint in the track fitting to improve the 
momentum resolution. This approach yields a very uniform tracking efficiency wifhin fhe accepfance, 
which is needed fo avoid geomefrical biases of fhe jef reconsfrucfion algorifhm caused by a non-uniform 
densify of reconsfrucfed fracks. The procedure is described in defail in fhe confexf of jef reconsfruc¬ 
fion wifh ALICE in Pb-Pb evenfs il2ll . Lor fhe analyzed dafa, fhe additional fracks (wifhouf SPD frack 
poinfs) consfifufe approximafely 4.3% of fhe used frack sample. Tracks wifh > 0.15 GeV/c and 
wifhin a pseudorapidify inferval |T]iab| < 0.9 are used as inpuf fo fhe jef reconsfrucfion. The overall effi¬ 
ciency for charged particle defecfion, including fhe effecf of fracking efficiency as well as fhe geomefrical 
accepfance, is 70% af px = 0.15 GeV/c and increases fo 85% af px = 1 GeV/c and above. 

2.2 Jet reconstruction and background corrections 


Lor the present analysis, the anti-kx algorithm from the LastJet package Il43ll has been used to reconstruct 
jets from measured tracks with resolution parameters of /? = 0.2 and R = 0.4. In general, jets are only 
considered for further analysis if the jet-axis is separated from the edge of the track acceptance in T]iab 
by at least the resolution parameter R used in the jet finding, e.g. jets reconstructed with R = 0.4 are 
accepted within |T]jet,iab| < 0.9 —0.4 = 0.5. The jet transverse momentum is calculated by LastJet using 
the px recombination scheme. To enable background corrections, the area A for each jet is determined 
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internally by distributing ghost particles into the area that is clustered t44l] . Ghost particles have vanish¬ 
ing momentum and therefore do not influence the jet finding procedure. By construction, the number of 
ghost particles in a jet is a direct measure for the jet area. A ghost particle density of 200 per unit area 
(0.005 area per ghost particle) was used to obtain a good area resolution with a reasonable computing 
time. 
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Fig. 1: (Color online) Left; Probability distribution of the event-by-event transverse momentum background den¬ 
sity (see Eq. [T]i. The mean and variance for two event classes are indicated in the figure. represents 

uncorrected jet pj. Right; Probability distribution of background fluctuations calculated with the random cone 
approach and defined via Eq. |4] (resolution parameter R = 0.4). 


In Pb-Pb collisions, the background from particles not from the same hard scattering as the jet has a 
significant impact on the reconstructed jet momentum 1121 13511. The transverse momentum density of this 
background is estimated with a statistically robust method by using the median of all jet px.chjet per area 
within one event for jets reconstructed with the kj algorithm. In p-Pb collisions, the multiplicity density 
is two orders of magnitude smaller than in central Pb-Pb collisions i40n . so a corresponding reduction 
of the jet background is expected. To obtain a reliable estimate for the more sparse environment of p-Pb 
events a modified version of the approach described in 1451] for pp collisions is employed. It uses the 
same method as in Pb-Pb, but contains an additional correction factor, C, to account for regions without 
particles, which otherwise would not contribute to the overall area estimate. The background density for 
each event is then given by 


Pch = median < 

^ I A,- 


•c, 


( 1 ) 


where i runs over all reconstructed kj jets in the event with momentum pTj and area A,. C is defined by 


C = 




( 2 ) 


Here, the numerator is the area of all jets containing tracks and the denominator, Aacc, is the acceptance 
in which charged particles are considered as input to the jet finding (2 x 0.9 x In). The probability 
distribution for the background density in this method, with the same track selection criteria as the signal 
jet reconstruction and a radius of 0.4, is shown in Fig. [T] (left). The background density obtained with 
R = 0.4 is used both for the correction of signal jets with R = 0.4 and R = 0.2 to avoid event-by-event 
fluctuations in the difference of the momenta for the two radii. 


The probability distribution of Pch decreases approximately exponentially. It is smaller than 4 GeV/c 
for 98.6% of all events. The mean background density and its variance for all events is (Pch) = 1-02 
GeV/c (with negligible statistical uncertainty) and cj(pch) = 0.91 ±0.01 GeV/c. For events containing 
a jet with uncorrected transverse momentum px.chjet > 20 GeV/c, it is (Pch) = 2.2 ±0.01 GeV/c and 
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(y{Pch) = E47 ±0.09 GeV/c, respeetively. The observed increase of the underlying event activity for 
events that contain a high-px jet is expected. This increase is already present in pp collisions and has 
been quantified in detail and with more differential observables than the background density, e.g. in 

The background density estimate provides an event-by-event correction for each jet with reconstructed 
transverse momentum px.chjet and jet area Ach jet: 


raw A 

PT,ch jet = Px.chjet ~^chjet ' Pch- 


(3) 


However, this approach neglects that the background for a given event is not uniformly distributed in 
the (Tjiab, <p)-plane but fluctuates from region to region. These fluctuations are mainly Poissonian, but 
also encode correlated region-to-region variations of the particle multiplicity and the mean px 13511 . The 
effect of these fluctuations can be accounted for on a statistical basis in the unfolding of the measured 
jet px,ch jet-distributions. The distribution of region-to-region density fluctuations around the event-wise 
background densiw estimate can be evaluated for the full event sample by a Random Cone (RC) approach 
as described in 113SH . Cones with a radius R corresponding to the resolution parameter of the jet finding 
algorithm are placed randomly in the (fliabj <P) jet-acceptance and the transverse momenta for all tracks 
(charged particles) falling into this cone are summed and compared to the background estimate: 


5pT,ch — ^PT.i — PchA, A — KR^. 


(4) 


The distribution of the residuals, 5px,ch> as shown in Fig. [T] (right) for R = 0.4, is a direct measure for 
all intra-event fluctuations of the background and can be used directly in the unfolding procedure. In 
Fig. [T] (right), a clear asymmetry of the distribution is visible. It is caused by the fact that the 5px,ch 
distribution of single particles sampled in the cone is asymmetric. Since the number of particles within 
a cone increases with its size, statistical fluctuations of the background estimate also increase (see also 
1 3j]). Furthermore, the randomly placed cones can also overlap with jets. In p-Pb collisions, there is the 
possibility for multiple hard collisions within one p-Pb event, so a jet can also be the background to a 
jet from another hard collision and contribute as an upward fluctuation. Therefore, an overlap of random 
cones with possible signal jets should not be a-priori excluded in the fluctuation estimate, but is part of 
its systematic uncertainty. 


2.3 Detector effects and unfolding 


The main detector-related effects on the reconstructed jet are the reconstruction efficiency and the mo¬ 
mentum resolution for single charged particles. To determine the correction for these, a full detector 
simulation of pp jet events generated with PYTHIA6 (Perugia 2011, version 6.425) |£7] and GEANT3 
particle transport ll^ is performed. In the simulation, two jet collections are matched geometrically 
(closeness in (T]iab, <p)-plane) with a one-to-one correspondence 11211 : jets reconstructed at the charged 
particle level (pari) without detector effects and jets reconstructed from tracks after particle transport 
through the AFICE detector (det). In the simulation, the particle level reconstruction includes charged 
primary particles produced in the collision with px > 0.15 GeV/c. Charged decay products from primary 
particle decays, excluding those from weak decays of strange particles, are included with the same px 
threshold. The response matrix is populated with matched particle- and detector-level jets. It relates the 
particle-level to the detector-level charged jet momentum and encodes the effects of single-particle mo¬ 
mentum resolution and reconstruction efficiency on the reconstructed jet momentum. A correction for 
the missing energy of neutral jet-constituents is not applied. The response is shown on a logarithmic scale 
in Fig. |2] (left) for charged jets with R = 0.4 and particle-level momentum between 45 < P^ljet < 50 
GeV/c. It can be seen that the most probable value for the reconstructed momentum is the particle-level 
momentum, but the distribution has large tails to the left and right. It is more probable that jets are 
reconstructed with a lower momentum than the truth, which is due to the dominating effect of the single¬ 
particle reconstruction efficiency that reduces the number of reconstructed particles in a jet. The tail to 
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Fig. 2: (Color online) Left; Projection of the combined unfolding matrix for jets with particle-level momentum 
45 < jet ^ GeV/c. The matrix is obtained from the combination of the detector response and background 
fluctuation matrices, which are also shown as projections (see text for details). Right: Probability distribution of 
the relative difference between particle-level (generated true) and detector-level charged jet transverse momentum 
for jets with different momenta. The effect of background fluctuations is included for the jets reconstructed at 
detector level. Characteristic values of the distributions are summarized in Table [T] 


the right-hand side is mainly due to the single-particle momentum resolution, where a fraction of tracks 
is reconstructed with higher momentum than the truth, causing an upward shift of the jet momentum. 


In addition, Fig.|2](left) shows the effect of the background fluctuations on the reconstructed jet momen¬ 
tum and the combination of detector effects and background fluctuations. Even though the background 
fluctuations show a strong tail to the right-hand side, it is seen that in the combined unfolding matrix 
the effects of single-particle momentum resolution play the dominant role in reconstructing a jet with 
momentum higher than the truth. The default algorithm for the unfolding of the measured jet spectrum 
is based on the Singular Value Decomposition (SVD) approach 11491] as implemented in the RooUnfold 
package 15011 . The default prior in the unfolding procedure is a smoothed version of the uncorrected jet 
spectrum itself. In addition to the SVD unfolding approach, Bayesian 115 iL 15211 and 15311 unfolding 
have been used for systematic comparisons and validity checks. The unfolded spectrum is also corrected 
for unmatched jets using a jet reconstruction efficiency obtained from generated-reconstructed compar¬ 
ison. This jet reconstruction efficiency is larger than 96% in the considered momentum range. 


The influence of these detector effects and background fluctuations on the jet momentum is shown for 
three transverse momentum intervals in Eig. |2] (right) via the probability distribution of the relative dif¬ 
ference of the detector-level and particle-level charged jet transverse momentum. Eor all momentum bins 
the distribution is asymmetric. The most probable response was determined using Gaussian fits to the 
peak region. It can be seen in Tab.[T]that it is close to zero (< 1%) with a mild px dependence. To further 
quantify the distributions, numerical values for their mean and width are also given in Tab. [T] Since the 
width is not a well-defined measure of the jet momentum resolution for these asymmetric distributions, 
the quartiles of the distribution are provided in addition. Approximately, 25% of the jets have a larger 
momentum than the generated. The 50% (median) correction is only 5% for = 20 — 25 GeV/c 

and increases towards larger jet momenta. In Tab. [T] the values for the respective distributions without 
background fluctuations are also given (not shown in Eig. |2ll. Clearly, the instrumental response domi¬ 
nates the jet response as already seen in Eig.[2](left). The main effect of the background fluctuations is a 
broadening of the jet response and an upward shift of the average reconstructed energy due to the asym¬ 
metric shape of the fluctuations as seen in Eig. [T] (right). The most probable value remains unaffected 
within the uncertainties when background fluctuations are included. 
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part 

Pj^ch jet 

20-25 GeV/c 

45 - 50 GeV/c 

80-90 GeV/c 

MPV (Gaussian Fit) 

0.006 ±0.002 

-0.001 ±0.002 

-0.010± 

0.004 

id. w/o bkg. fluct. 

0.007 ±0.001 

-0.003 ±0.002 

-0.013 ± 

0.004 

Mean 

0.149 ±0.030 

-0.181 ±0.030 

-0.222 ± 

0.030 

id. w/o bkg. fluct. 

-0.163 ±0.030 

-0.188±0.030 

-0.226 ± 

0.030 

Width a 

0.238 ±0.030 

0.246 ±0.030 

0.259± 

0.030 

id. w/o bkg. fluct. 

0.233 ±0.009 

0.245 ±0.005 

0.258 ± 

0.003 

Quartile, 25% above 

0.01 ±0.01 

-0.01 ±0.02 

-0.01 ± 

0.02 

id. w/o bkg. fluct. 

0.01 ±0.01 

-0.01 ±0.02 

-0.03 ± 

0.02 

Quartile, 50% above 

-0.05 ±0.04 

-0.09±0.01 

-0.13± 

0.04 

id. w/o bkg. fluct. 

-0.07 ±0.04 

-0.09 ±0.04 

-0.13± 

0.04 

Quartile, 75% above 

-0.25 ±0.06 

-0.29 ±0.05 

-0.37 ± 

0.04 

id. w/o bkg. fluct. 

-0.27 ±0.04 

-0.29 ±0.06 

-0.37 ± 

0.04 


Table 1: Characteristic values for the distribution of residuals of the total charged jet response shown in Fig. |2] 
(right), including the effect of background fluctuations and without: most probable value (MPV) determined via a 
Gaussian fit to the central peak region, first and second moment (mean and width a), and quartiles. The precision 
of the quartiles is limited by the finite bin width of 0.01. 

2.4 Nuclear modification factor 

The nuclear modification factor compares a pj-differential yield in p-Pb collisions to the differential 
production cross section in pp collisions at the same -^/^nn to quantify nuclear effects: 

_ d^ApPb/dT7dpT 
(EpPb) • d2 app /dT] dpT' 


Here, (TpPb) is the nuclear overlap function which accounts for the increased parton flux in p-Pb com¬ 
pared to pp collisions. It is related to the number of binary nucleon-nucleon collisions via (TpPb) = 


(Ncoii) /tt'jNEL ^'^d has been calculated in a Glauber Monte Carlo, as described in ll38l] . Here, rep¬ 


resents the total inelastic cross section in pp collisions. For minimum bias p-Pb collisions, the nuclear 
overlap function is (TpPb) = (0.0983 ±0.0034) mband (Ncoii) = 6.87 ±0.56. In this paper, the refer¬ 
ence differential production cross section in pp is constructed from the ALICE charged jet measurement 
at 7 TeV by a pQCD based scaling. In the nuclear modification factor, the invariant yield for NSD 
events in p-Pb is compared to inelastic pp collisions. Hence, the additional correction of (2.3 ± 3.1)% is 
applied as discussed above. 

2.5 NLO calculations and pp reference 

Perturbative QCD calculations are used for two purposes in this paper: for comparison to the measure¬ 
ment of jet production in p-Pb, and as additional input to the construction of the pp reference. The 
calculations have been performed within the POWHEG box framework iSSl. 15611 . which facilitates next- 
to-leading order (NLO) precision in calculating parton scattering cross sections in an event-by-event 
Monte Carlo. Event-by-event the outgoing partons from POWHEG are passed to PYTHIA8 i57n where 
the subsequent parton shower is handled. For this, a POWHEG version matched to the PYTHIA8 frag¬ 
mentation is used to avoid double counting of NLO effects already considered in the PYTHIA8 code. 
The Monte Carlo approach has the advantage that the same selection criteria and jet finding algorithm 
can be used on final state particle level, as in the analysis of the real data, in particular, the limitation to 
charged constituents of a jet. The dominant uncertainty in the parton level calculation is given by the 
choice of renormalization scale, /Tr, and factorization scale, /Tp- The default value has been chosen to 
be /Tr = /^F = Pt and independent variations by a factor of two around the central value are considered 
as the systematic uncertainty. In addition, the uncertainty on the parton distribution functions has been 
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taken into account by the variation of the final results for the respective error sets of the parton density 
functions (PDFs). 


For the comparison with the measured p-Pb data, proton PDFs corrected for nuclear effects (CTEQ6.6 
1 581] with EPS09 115911 ) have been used. Prior to passing the scattered partons to PYTHIAS for showering, 
they can be boosted into the same reference frame as the p-Pb reaction by yNN = 0.465. 


The construction of the pp reference at = 5.02 TeV is based on the AFICE measurement of charged 
jets in pp collisions at 7 TeV, described in detail in 15411 . For the purpose of the reference scaling, the 
same analysis chain has been used as for p-Pb. The same binning in pseudorapidity and transverse 
momentum allows for a partial cancellation of common systematic uncertainties in the pp and p-Pb data 
sets. In addition, the same background subtraction approach as in the p-Pb analysis is used for the pp 
data. In the present analysis, the scaling is done with a factor which is determined for each px.chjet bin 
by the NFO pQCD calculations (POWHEG+PYTHIA8) at the two energies. For the pp reference scaling 
the parton distribution functions in the POWHEG calculation have been replaced by the free proton PDF 
from CTEQ6.6. The scaling factor is given by 


P{Pi) = 


yield(pT,chjet)| 


5.02 TeV,boosted 
pp,NLO 


1 7 TeV 


yield(PT,chjet)|pp_NLO 


( 6 ) 


The factor decreases monotonically from F k, 0.65 to 0.45 in the reported pj range. As already described 
above, the laboratory frame is not the center-of-mass frame of the collision as is the case for pp collisions. 
Therefore, the numerator of the scaling factor F in Eq.^is determined in the NFO calculation where the 
additional Forentz-boost is applied to the hard scattered partons prior to fragmentation. The resulting 
reduction of the observed jets for |T]iab| < 0.5 is smaller than 5% in the relevant momentum range. 


2.6 Jet production cross section ratio 


The broadening or narrowing of the parton shower with respect to the original parton direction can have 
a direct impact on the jet production cross section reconstructed with different resolution parameters. 
This can be tested via the ratio of yields or cross sections in common rapidity interval, here Ifiiabl < 0.5 
for R = 0.2 and 0.4: 


.^(0.2, 0.4) 


dCJpPb,R=0.2/dpT 

dCJpPb,R=0.4/dpT 


(V) 


Considering the extreme scenario that all fragments are already contained within /? = 0.2 this ratio is 
unity. In this case, also the statistical uncertainties between R = 0.2 and R = 0.4 are fully correlated and 
cancel completely in the ratio, when the jets are reconstructed from the same data set. In the case the 
jets are less collimated, the ratio decreases and the statistical uncertainties only cancel partially. For the 
analysis presented in this paper, the conditional probability for reconstructing a 7? = 0.2 jet in the same 
px-bin as a geometrically close R = 0.4 jet is 25 — 50%, which leads to a reduction of the statistical 
uncertainty of the ratio of 5 — 10% compared to the case of no correlation. 


2.7 Systematic uncertainties 

The various sources of systematic uncertainties are listed in Tab. |2] for the full px-range of the three 
observables presented in this paper: jet production cross section, nuclear modification factor, and cross 
section ratio. The most important sources will be discussed in the following. 

The dominant source of uncertainty for the px-differential jet production cross section is the imperfect 
knowledge of the single-particle tracking efficiency that has a direct impact on the correction of the jet 
momentum in the unfolding, as discussed above. In p-Pb collisions, the single-particle efficiency is 
known with a relative accuracy of 4%, which is equivalent to a 4% uncertainty on the jet momentum 
scale. To estimate the effect of the tracking efficiency uncertainty on the jet yield, the tracking efficiency 
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is artificially lowered by randomly discarding a certain fraction (4% in p-Pb) of tracks used as input for 
the jet finding. Depending on the shape of the spectrum, the uncertainty on the single particle efficiency 
(jet momentum scale) translates into an uncertainty of 8 to 15% on the yield. 

To estimate the uncertainty on the p-Pb nuclear modification factor, the uncertainty on the single-particle 
tracking efficiency in the two collision systems (pp and p-Pb) has to be evaluated. This uncertainty on the 
efficiency is correlated between the data sets, since the correction is determined with the same underlying 
Monte Carlo description of the ALICE detector and for similar track quality cuts. Only variations of de¬ 
tector conditions between run periods may reduce the degree of correlation. The uncorrelated uncertainty 
has been estimated to be 2%, and the uncertainty for the nuclear modification factor has been determined 
by artificially introducing such a difference in the tracking efficiency between the two collision systems. 

The uncertainty on the spectra induced by the underlying event subtraction has been estimated by com¬ 
paring the results with various methods for background subtraction; ranging from purely track-based to 
jet-based density estimates, including an T]iab-dependent correction. As seen in Eig. [T] a typical correc¬ 
tion nR^Pch for a jet with R = 0.4 is about 1 GeV/c. The uncertainty on this correction can be treated 
similar to an uncertainty on the jet momentum scale. Eor the final spectrum, the uncertainty on the yields 
from the background correction method is approximately 2%. 

In the determination of the fluctuations of the underlying event, the main uncertainty is given by the 
exclusion of reconstructed jets in the random cone sampling of the event. The probability for a random 
cone to overlap with reconstructed jets is higher than for the jets itself. On average, a jet can overlap with 
^coii — 1 jets in one event. The random cone can overlap with Acoii jets. To account for this, the Spx.ch 
calculation can be modified to discard on a statistical basis random cones that overlap with signal jets. 
This lowers the average overlap probability. However, since this modified 5pT,ch calculation strongly 
depends on the signal jet definition and also on how an overlap is defined, it is not used by default but 
considered for systematic uncertainties. The effect of this partial signal exclusion approach on the fully 
corrected jet yields is of the order of 0.1%. 


The uncertainty of the scaling procedure to obtain the reference spectrum is estimated by determining 
the scaling factors E(pt) after varying the scales Hr and /Tp in the POWHEG NLO generation, and by 
using different tunes in the outgoing fragmentation handled by PYTHIAS. Eurthermore, standalone cal¬ 
culations with PYTHIA6 and PYTHIAS using different generator tunes and with HERWIG at the two 
energies have been performed to obtain scaling factors according to Eq. 0 A general uncertainty for 
how well LO generators and NLO calculations can describe the v^-dependence of particle production 
is also considered: in ALICE measurements of the production in pp collisions, it has been observed 
that pQCD calculations predict a stronger increase of the production cross section when going from 0.9 
to 7 TeV than supported by the data ifiOtl . A similar effect is also seen in unidentified charged hadrons 
measured with ALICE at 0.9, 2.76, and 7 TeV 16IH . Eurthermore, the -y^-dependence of the jet produc¬ 
tion cross section has been cross checked internally with an interpolation between 7 and 2.76 TeV, using 
preliminary ALICE results on charged jets at = 2.76 TeV. In total, these studies yield an additional 
uncertainty on the pp reference of 10% for the extrapolation from 7 to 5.02 TeV. It is reported as an 
independent normalization uncertainty, similar to the uncertainty on the nuclear overlap function. 


3 Results 

The pt- differential production cross sections for jets reconstructed from charged particles in mini¬ 
mum bias p-Pb collisions at y/^NN = 5.02 TeV are shown in Pigs. [3] and |4] for the resolution param¬ 
eters R = 0.4 and R = 0.2. The spectra are found to agree well with scaled NLO pQCD calculations 
(POWHEG-i-PYTHIAS) using nuclear PDPs (CTEQ6.6-I-EPS09) as seen best in the ratio data over calcu¬ 
lation in the lower panels. However, the effect of the nuclear PDPs on the jet production in the reported 
kinematic regime is almost negligible, as seen in the comparison to calculations with only proton PDPs 
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Fig. 3: (Color online) Top panel; px-differential production cross section of charged jet production in p-Pb colli¬ 
sions at 5.02 TeV for R = 0.4. Bottom panel; Ratio of data and NLO pQCD calculations. The global uncertainty 
from the measurement of the visible cross section of 3.5% is not shown. The uncertainties on the pQCD calculation 
are only shown in the ratio plot as dashed lines. The pQCD calculations take into account the rapidity shift of the 
nucleon-nucleon center-of-mass system in p-Pb with a boosted parton system. 
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Fig. 4: (Color online) Top panel; px-differential production cross section of charged jet production in p-Pb colli¬ 
sions at 5.02 TeV for R = 0.2. Bottom panel; Ratio of data and NLO pQCD calculations. The global uncertainty 
from the measurement of the visible cross section of 3.5% is not shown. The uncertainties on the pQCD calculation 
are only shown in the ratio plot as dashed lines. The pQCD calculations take into account the rapidity shift of the 
nucleon-nucleon center-of-mass system in p-Pb with a boosted parton system. 
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Observable 

Jet cross section 

R 

pPb 


Resolution parameter 

Uncertainty source 

R = 0.2 

R = 0.4 

R = 0.2 

R = 0A 

0.2/0.4 

Single-particle efficiency (%) 

7.9-12.8 

10.2-14.2 

4.1-5.9 

4.9-6.3 

2.1-2.1 

Unfolding (%) 

2.2 

1.7 

2.8 

2.2 

1.5 

Unfolding prior steepness (%) 

1.4-4.8 

0.5-4.0 

2.9-8.0 

0.9-4.4 

1.1-1.5 

Regularization strength (%) 

3.1-3.9 

2.3-4.4 

CO 

1 

2.3-5.6 

1.1-4.7 

Minimum px cut-off (%) 

1.1-0.3 

2.3-0.1 

1.3-1.4 

2.8-4.1 

1.2-0.4 

Background estimate (%) 

1.8-0.6 

3.7-1.5 

1.8-0.6 

3.7-1.5 

2.0-0.9 

5px,ch estimate (%) 

0.0-0.0 

0.1-0.0 

0.0-0.0 

0.1-0.0 

0.1-0.0 

Combined uncertainty (%) 

9.2-14.4 

11.5-15.5 

7.1-11.9 

7.5-10.7 

3.8-5.7 

(Vb) (%) 

- 

- 

3.4 

3.4 

- 

pp cross section (%) 

- 

- 

3.5 

3.5 

- 

Reference scaling pp 7 TeV (%) 

- 

- 

10.0 

10.0 

- 

NSD selection efficiency p-Pb (%) 

- 

- 

3.1 

3.1 

- 

Combined scaling uncertainty (%) 

- 

- 

11.6 

11.6 

- 


Table 2: Summary of systematic uncertainties on the fully corrected jet spectrum, the corresponding nuclear 
modification factor, and the jet production cross section ratio for the resolution parameters R = 0.2 and R — 0.4. 
The percentages are given for the whole shown transverse momentum range 20 — 120 GeV/c 


(CTEQ6.6). 

Eigure m also shows the jet spectra for —0.65 < Tjiab < —0.25 and 0.25 < T]iab < 0.65 compared to the 
results from the symmetric selection |T]iab| < 0.5. Here, T]iab denotes the pseudorapidity of the jet axis. 
The first selection roughly corresponds to a small window around mid-rapidity for the nucleon-nucleon 
center-of-mass system, while the second is separated from it by about one unit in rapidity. No significant 
change of the jet spectra is observed for these two T]iab regions centered at —0.45 and 0.45. Thus, the jet 
measurement has no strong sensitivity to the rapidity shift and the pseudorapidity dependent variation of 
the multiplicity (underlying event) within the statistical and systematic uncertainties of the measurement. 


The nuclear modification factor /?pPb is constructed based on the px-differential yields and the extrapo¬ 
lated pp production cross section at 5.02 TeV for R = 0.2 and 0.4. It is shown in the left and right panel 
of Eig. m respectively. In the reported px-range, it is consistent with unity, indicating the absence of a 
large modification of the initial parton distributions or a strong final state effect on jet production. Be¬ 
fore comparing these results to the measured single-particle results for/?pPb, one has to consider that the 
same reconstructed px corresponds to a different underlying parton transverse momentum. Assuming 
that all spectra should obey the same power law behavior at high px, an effective conversion between the 
spectra can be derived at a given energy via the POWHEG-I-PYTHIA8 simulations described above. To 
match the single charged particle spectra in the simulation to charged jets with R = 0.4, a transformation 
^ 2.2Sptf is needed. Thus, the reported nuclear modification factor for charged jets probes roughly 
the same parton px-region as the AEICE measurement of single charged particles that shows a nuclear 
modification factor in agreement with unity in the measured high-px range up to 50 GeV/c 12711 . 


Since the jet measurements integrate the final state particles, they have a smaller sensitivity to the frag¬ 
mentation pattern of partons than single particles. Differences between the nuclear modification factor 
for jets and single high-px particles, as suggested by measurements in |28, 2^, could point to a modified 
fragmentation pattern or differently biased jet selection in p-Pb collisions. 


A modified fragmentation pattern may be also reflected in the collimation or transverse structure of jets. 
A first step in testing possible cold nuclear matter effects on the jet structure is the ratio of jet production 
cross sections for two different resolution parameters. It is shown for R = 0.2 and R = 0.4 in p-Pb in 
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Fig. 5: (Color online) Nuclear modification factors Rpp\, of charged jets for R = 0.2 (left) and R = 0.4 (right). 
The combined global normalization uncertainty from (TpPb), the correction to NSD events, the measured pp cross 
section, and the reference scaling is depicted by the box around unity. 
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Fig. 6: (Color online) Charged jet production cross section ratio for different resolution parameters as defined 
in Eq. [T] The data in p-Pb collisions at i/snn = 5.02 TeV are compared to PYTH1A6 (tune; Perugia 2011, 
no uncertainties shown) and POWHEG+PYTHIA8 (combined stat. and syst. uncertainties shown) at the same 
energy, and to pp collisions at 7 TeV (only stat. uncertainties shown). 
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Fig. 0 and compared to PYTHIA6 (Tune Perugia 2011) and POWHEG+PYTHIA8 at y^nn = 5.02 TeV 
and to AEICE results in pp collisions at ^/s = 7 TeV 15411 . All data show the expected increase of the 
ratio from the increasing collimation of jets for higher transverse momentum and agree well within the 
uncertainties. No significant energy dependence or change with collision species is observed. The data 
for p-Pb collisions is well described by the NEO calculation as well as by the simulation of pp collisions 
with PYTHIA6 at the same energy. It should be noted that the ratio for charged jets is, in general, above 
the ratio obtained for fully reconstructed jets, containing charged and neutral constituents. This can be 
understood from the contribution from neutral pions that decay already at the collision vertex and lead to 
an effective broadening of the jet profile when including the neutral component in the jet reconstruction, 
mainly in the form of decay photons. For the same reason, the inclusion of the hadronization in the NEO 
pQCD calculation is essential to describe the ratio of jet production cross section as also discussed in 


4 Summary 

In this paper, px-differential charged jet production cross sections in p-Pb collisions at = 5.02 

TeV have been shown up to px.chjet of 120 GeV/c for resolution parameters R = 0.2 and R = 0.4. The 
charged jet production is found to be compatible with scaled pQCD calculations at the same energy using 
nuclear PDFs. At the same time, the nuclear modification factor 7?pPb (using a scaled measurement of jets 
in pp collisions at ^/s = 7 TeV as a reference) does not show strong nuclear effects on jet production and 
is consistent with unity for R = 0.4 and R = 0.2 in the measured px-range between 20 and 120 GeV/c. 
The jet cross section ratio of /? = 0.2/0.4 is compatible with 7 TeV pp data and also with the predictions 
from PYTHIA6 Perugia 2011 and POWHEG+PYTHIA8 calculations at 5.02 TeV. No indication of a 
strong nuclear modification of the jet radial profile is observed, comparing jets with different resolution 
parameters R = 0.2 and R = 0.4. 
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